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1. Introduction 

Dynamical overlap simulations are extremely time consuming [1]. First large scale results 
have been reported by JLQCD at this conference, but very large computer resources were needed 
to obtain them. An efficient alternative to fully dynamical simulations with chiral fermions might 
be provided by the so-called mixed action approach. Here, the fermion discretization employed for 
the sea quarks is different from the one used for the valence quarks. Provided the sea and valence 
quark masses are properly matched such a setup is expected to give the same continuum physics as 
standard unquenched lattice QCD. 

The main idea is to use sea quarks which are (relatively) cheap to simulate and therefore large 
volumes, good statistics and small quark masses are attainable even with the fermion determinant 
included. Such fermion discretizations, however, typically break chiral symmetry which makes 
the extraction of certain observables difficult. These difficulties can be reduced by using valence 
quarks that fulfill the Ginsparg-Wilson equation [2] and thus respect chiral symmetry at finite 
lattice spacing, e.g. Neuberger's overlap fermions [3], domain wall fermions [4, 5] with a large 
extent of the fifth dimension or perfect actions [6]. 

One drawback of the mixed action approach is that the resulting theory is not unitary at non- 
zero lattice spacing. This introduces well known diseases into the theory which, however, are 
expected to vanish in the continuum limit. Moreover, the effects of unitarity violation can be 
described by mixed action chiral perturbation theory [7, 8]. 

At this conference, we present results obtained on configurations generated by the European 
Twisted Mass Collaboration (ETMC) [9, 10]. We use an ensemble with two flavors of twisted mass 
fermions on a lattice of size 16^ x 32 and a lattice spacing of about a 0.1 fm. The (charged) 
pseudo-scalar mass is about mps ~ 300 MeV. 

2. Setup 

The lattices are generated with a tree-level Symanzik improved gauge action at j3 = 3.9. Ex- 
tensive studies have shown that this gauge action is a good compromise between ease of simulation 
and good short distance behavior [11]. 

2.1 Twisted mass sea fermions 

We take two twisted mass sea fermions into account. The corresponding Dirac operator reads 

£>tm = D^w^^^^ (mo) + ms T3 (2. 1) 

with D^/v the standard Wilson Dirac operator taken at un-twisted mass parameter mq. For more 
details see Refs. [12, 13]. In our simulation we set K = 0.1609 which is close to but not quite 
the critical value of k" = 0.160856. The twisted mass parameter (which determines the bare quark 
mass) is set to ajj. = 0.004. 

Twisted mass QCD has many appealing properties. The fermion measure is finite and positive 
everywhere and, if when K is set to its critical value, physical observables are automatically ff{a) 
improved [13]. ETMC has demonstrated that twisted mass fermions can be simulated at small 
pseudo-scalar masses and large volume with currently available computer resources, see the talks 
by Jansen and Urbach at this conference [9] and Shindler at ICHEP [10]. 
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Figure 1: The effect of HYP smearing on the low-lying eigenvalues of the kernel operator /!(—7?o)- We plot 
the average of the magnitude of 20th eigenvalue as a function of the negative mass shift — ^o- This value 
gives the lower bound of the rational approximation used in the construction of the overlap. The data is from 
an ensemble at j3 = 3.75, k = 0.1661 and api = 0.005. 



2.2 Overlap valence fermions 

We use overlap fermions in the valence sector. They have exact chiral symmetry at finite 
lattice spacing. This ensures ^(a) improvement at the level of the action and of the suitably chosen 
operators. Since also the sea quarks are &{a) improved, we expect small scaling violations for the 
whole set-up too. The overlap operator is constructed from the Hermitian Wilson Dirac operator 
h = ysD\f which is taken at a negative mass shift —Rq. 

D,,im) = {Ro - y) [1 + 75e{h{-Ro))] +m (2.2) 

with £ the matrix sign function and m the bare quark mass. The parameter Rq is a mass at the cut- 
off scale which has to be tuned. The disadvantage of these fermions is the cost of the construction 
of the sign function. This turns out to be ^(100) times more expensive than the application of the 
kernel operator. 

We use the standard Hermitian Wilson Dirac operator h{—Ro) as kernel for the overlap, which 
we construct from HYP blocked links [14]. Among the benefits of the HYP blocking are lower 
cost of construction and better locality of the overlap operator (see Fig. [l]). We investigated one 
and two levels of stout smearing. In particular two steps can lead to a comparable situation for the 
overlap operator. However, the smaller range and greater (generally positive) experience with the 
HYP blocking lead us to use it in our production runs. 

It now remains to determine the best value of the negative mass shift —Rq of the kernel opera- 
tor. We set it such that the resulting overlap operator has optimal locality. Locality means that the 
the norm of the operator has an exponential decay at large distances 

\\D,,{x,y;U)\\<Ce-''^'''y^ . 



3 



Overlap fermions on a twisted mass sea 



S. Schaefer 




Figure 2: The locality test as defined in Eq. 2.3 We observe an exponential decay with a decay rate 
V ~ 0.607 for the optimal choice /?o = 1- 



To verify this, we follow Ref. [15] and measure the exponential decay of f{r) which is essen- 
tially the norm of overlap operator on a delta source. 



\ir = Do^r\ with t] (x) = ^^-.^o 
f{r) =max{||vA(x)|| : dist(x,xo) 



r] 



(2.3) 



with dist(x,y) the taxi-driver distance between x and y. The result is shown in Fig. ^for various 
choices of the parameter Rq. We find that the overlap operator has best locality properties if we set 
Rq = \, the free field value. (This is a well known effect of the smoothing of the gauge fields by 
HYP blocking.) For this choice of Rq we measure a decay rate v = 0.607. 

2.3 Matching 

The last parameter to fix is the valence quark mass. A simple matching condition is given by 
setting the masses of the pseudo-scalar mesons made of of two sea and two valence quarks equal, 



M 



ss 



vv ■ 



(2.4) 



For twisted mass fermions, the neutral and the charged pseudo-scalar are not degenerate. We chose 
to match the charged meson mass because of its smaller statistical uncertainty. In Fig. ^ we show 
the valence pseudo-scalar mass as a function of the bare quark mass. This has to match the meson 
mass from the sea quarks whose la band is indicated by the dashed lines. The matching quark 
mass parameter is close to aniq = 0.01. 



3. WChPT for twisted mass LQCD 

Mixed action QCD can be described by mixed ChPT provided the pseudoscalars are suffi- 
ciently light [7]. Mixed ChPT is formulated analogously to ChPT for lattice theories with the same 
type of sea and valence quarks. The chiral Lagrangian through O(a^) for untwisted Wilson sea 
and GW valence quarks, which is easily generalized to twisted Wilson sea quarks, can be found in 
Ref. [16]. 
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Figure 3: Matching the mass of the charged pseudo-scalar meson from the valence and the sea quarks. The 
dashed lines indicate the 1(7 range from the sea sector. 



At leading order one finds the following expressions for the mesons masses (my , jXs denote 
renormalized quark masses, assuming maximal twist in the sea sector). 

Val-Val: M^y = IBmy 

Sea-Sea: Mf^ = 2B^s, ^ig > 



Val-Sea: M^^ = B{mv + Ms) + a^Q 



Mix 

The valence-valence pion mass vanishes for my = because of the exact chiral symmetry in the 
valence sector. The minimal sea-sea pion mass mJ^ = 2\c2\a^ is proportional to the low-energy 
constant C2 [17] which determines the phase diagram for the Wilson sea quarks [18, 19]. Note that 

because of the extra low-energy constant Cmix present in mixed theories, and a measurement of the 
mixed meson mass is a direct measure for the size of the O(a^) cut-off effects in the mixed theory. 

Since mixed action theories are not unitary at non-zero lattice spacing, they are expected to 
show the same pathologies as partially quenched theories. A sensitive probe for these pathologies 
is the flavor non-singlet scalar correlator 

C{t) = Y,Wu{x,t)ud(Si,0)\<d) . 

X 

Using the spectral decomposition the scalar correlator can be written as 

C{t) = Ae-^-o' +B{t) + ..., (3.1) 

where Ma^ denotes the mass of the ao, the lightest 7=1 scalar meson. B{t) denotes the contribution 
of two pseudoscalar states (e.g. ttt]'), while the ellipsis represents excited scalar states, multi- 
hadron states, etc. which are thought to be less important and therefore ignored. 
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Figure 4: The scalar correlator for four values of the valence quark mass. The dotted lines represent the 
chiral perturbation theory prediction with the parameters determined from a fit to the data. 



The contribution B{t) from two pseudoscalar states has been computed in leading order (LO) 
ChPT including the leading cut-off effects of O(a^) [20, 21]. For large times t one finds approxi- 
mately {Nf = l, Mss = Mvvf 



B{t) 



1 + ^. 

Mvv 



(3.2) 



The term proportional to Jssa^t dominates the large t behavior. Note that, depending on the sign 
of the low-energy constant 755, the scalar correlator can become negative even if the matching 
condition M55 = Myy is chosen, in contrast to continuum partially quenched ChPT. 

Figure ^ shows the scalar correlator for various valence quark masses, together with a com- 
bined fit to the ChPT prediction.^ The scalar correlator is indeed negative for all valence quark 
masses shown. Even though ChPT is able to qualitatively describe the negative correlator, the fit is 
poor. One reason might be that the NLO correction are not negligible. 



4. Summary 

We presented first results from a mixed action project where we use the overlap Dirac oper- 
ator in the valence sector and twisted mass fermions in the sea. We discussed the choices of the 
parameters of the overlap operator: the smearing of the gauge links for expedience and improved 
locality of the operator, the tuning of the radius of the Ginsparg-Wilson circle Rq for optimal lo- 
cality and the determination of the valence quark mass such that the pion mass in valence and sea 
sector match. 

As physics result we discussed the scalar correlator and how to describe it with mixed action 
chiral perturbation theory. We found good qualitative agreement. 



' For simplicity we have sent the singlet mass parameter mo to infinity, see Ref . [20] for its definition. 
^For the fit we use the exact formula for C{t) with the singlet parameter mo = IGeV, not eq. (3.2). 
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